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Information on the intrinsic deformation of the proton can be obtained by studying 
the electromagnetic p — » A+ quadrupole transition. Recent experiments have 
shown that the electric quadrupole {E2) strength in 7p — > A+ is about 10 times 
larger than predicted by the simple quark model using only one-body currents. 
Our analysis provides evidence for the dominance of exchange currents in the 
Af — > A quadrupole transition, and identifies the physical mechanism leading to 
the observed E2 strength. 



1 Introduction 

The quadrupole moment of a particle measures the deviation of its internal 
charge distribution from spherical symmetry. However, as a particle with 
total angular momentum J = 1/2, the nucleon does not have a quadrupole 
moment in the laboratory frame. In order to learn something about the shape 
of the nucleon one has to electromagnctically excite it, e.g. to the A(1232) 
resonance, with total angular momentum J — 3/2, or to higher resonances. 

There are three different ways to electromagnctically produce a A(1232). 
Aside from the dominant ^ A magnetic dipole (Ml) excitation mode, in 
which the spin of a single quark is flipped, transverse electric quadrupole (£'2), 
and longitudinal charge quadrupole (C2) transitions are allowed by angular 
momentum and parity selection rules. The strengths of these electromag- 
netic multipoles at photon three-momentum transfer q = 0, called Gmi(0), 
G_B2(0), and Gc2(0), can be extracted from high precision photo-pion- and 
electro-pionproduction experiments off the proton H. The quadrupole transi- 
tion amplitudes Ge2{0) and Gc2(0) are a measure of the intrinsic defomiation 
of the nucleon. The empirical E2 strength is with G^l* — 0.133(20) EJ about 
10 times larger than predict/^f^y the first quark model calculations using 
spatial single quark currents 
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2 Single Quark — > A Quadrupole Transition 



Early quark model calculations of the ^ A quadrupole transition used 
aside from a spin-independent confin£Lnient a one-gluon exchange potential 
yOGEP i^g^-^ggn constituent quarks Ha. Without retardation and for equal 
quark masses it reads 
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where r = r j — r^ ; ai is the Pauli spin matrix, and is the color operator of 
the i-th quark. Usually, a regularized one-pion exchange interaction between 
constituent quarks is introduced, in order to account for the chiral symmetry 
of QCD and its spontaneous breaking at low energies 
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Here, fi is the pion mass, A the cut-off mass, and Ti denotes the isospin of the 
i-th quark. 

The tensor terms in the effective gluon and pion exchange potentials in- 
duce D-wave admixtures in the three-quark N and A wave functions. For 
clarity we restrict ourselves to a two-state harmonic oscillator model for the 
wave functions: 

\N) ^ as \{S ^ 1/2, 0) J = 1/2) + \{S = 3/2 , L = 2) J = 1/2) 
|A) = 6s 1(5 = 3/2 , L - 0) J = 3/2) + boKS = 1/2 , L = 2) J = 3/2) , (3) 

where the Z?-states are of mixed symmetric type. In Eq.(|3|), the inner spin S 
couples with the orbital angular momentum L to the total angular momentum 
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Table 1. Admixture coefficients for the S and mixed symmetric D states in the nucleon 
and A ground state wave function as calculated by different authors. 
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0.056 


Ref. 






0.934 


-0.047 
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0.064 



J of the auarks. The D-wave probabihties aj^ and bj^ as calculated by different 
authors aa lie between 0.2% and 0.4% (see table |]). Ref.H obtains a larger In- 
state probability of ~ 1%. In any case, the calculated D state probabilities are 
much smaller than the 5% Z?-state admixture in the deuteron wave function. 
We will come back to this point. 

In the single quark transition model, based on the spatial one-body cur- 
rent J[x], a quark moving in an excited D state in the nucleon can absorb 
electromagnetic quadrupolc radiation (L-y—2), and thus fall to the S state in 
the A. The amplitude for this process is proportional to aubs- In addition, a 
quark in an S state in the nucleon can absorb quadrupole radiation (L^=2), 
and jump to an excited D-state in the A. This happens with an amplitude 
asbjj. Both transitions, graphically displayed in Fig. change the orbital an- 
gular momentum L, but leave the inner spin state S of the quarks unchanged. 
The resulting single-quark E2 transition strength is E 



G_E2(J[i]) = {aobs + asbn) ■ 



(4) 



The Ge2 used here is dimensionless. The connection with the dimensionful 
transition quadrupole moment, or the helicity amplitudes is given in Ref. 13. 
With the D state admixtures calculated with the iffiuripus potentials based on 
Eq.(|l]) and Eq.(||) one obtains Ge2 = 0.003 — 0.01 BHla, i.e. theoretical results 
that differ by an order of magnitude from the recent experimental results. 
This suggests that some important dynamical feature is still missing. 
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Figure 1. Conventiptial explanation of the 7A'^ A quadrupole transition via the one-body 
current J[x] in FigMa) (impulse approximation). In this approximation the E2 transition 
amplitude is a coherent superposition of the two L changing but S conserving one-body 
transitions (upper and lower part of the figure). In a single-quark transition, the absorption 
of a E2 photon is therefore only possible if either the nucleon (left) or the A (right) contains 
a D wave admixture (deformed valence quark orbit). The transition probability amplitude 
GE2{-f[i]) of Eq.(W) is strongly suppressed due to the small D wave admixtures in the N 
und A wave function. 



3 Exchange Currents and Siegert's Theorem 

From the theoretical point of view the conventional explanation outlined above 
is incomplete because it violates the continuity equation: 

q-J{q) = [H,p{q)]. (5) 

After a decomposition of the charge and current operators into one- and two- 
body terms, i.e., p = p^^ + /9[2], and J = J[i] -|- J[2], one can show that the 
spatial two-body current Jp] satisfies the following consistency relation 

q- J[2](q) = [V[2], P[i](q)] • (6) 

Equation (^) connects the quark-quark interactions V[2] , which determine the 
coefficients as, an, bs, bo in Eq. (^) with the two-body currents J[2] of Fig.(||), 
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Figure 2. Feynman diagrams of the four vector current J** = (p, J): (a) one-body current 
Jj'j^j , and (b-d) two-body gluon and pion exchange currents J . If the quarks interact via 
pion and gluon exchange, the two-body exchange currents depicted in diagrams (b-d) must 
be taken into account. If only diagram (a) is considered, the continuity equation is 
violated. 



which determine the electromagnetic properties of the N—A system. Eq.(|6|) is 
violated if the potential contains momentum and/or isospin dependent terms 
but the electromagnetic current contains only one-body terms. In fact, the 
potentials of Eq.(|l|) and Eq.(||) do not commute with the one-body charge 
operator, thus implying a two-body current J^j- m 

An important theorem based on Eq.(|^) is Siegert's theorem □. In the 
limit of small momentum transfers it relates the transverse electric and 
longitudinal Coulomb multipoles T*-^'^: 

(/|T^^(|q| - 0)1*) = (/|T^^(|q| - 0)^) . (7) 

Thus, Ge2 can be calculated via the charge operator (right hand side). The 
result will be the same as the one based on the spatial current operator (left- 
hand side). 

However, until recently, calculations of the E2 transition strength Ge2 
using the spatial one-body current operator have differed considerably (in 
some calculations by an order of magnitude) from those using the one-body 
charge jDperator and Siegert's theorem. This important observation was made 
in Ref.tl, and has been confirmed by other authors Several improvements 
have been proposed, e.g. an increase of the number of harmonic oscillator 
states, and the inclusion of relativistic corrections to the one-body current, in 
order to remove this difference. However, the main reason for this discrepancy 
has not been recognized in previous works. We have recently shown Q that this 
difference is almost entirely explained by spatial two-body exchange currents 
J [2] required by Eq.(^) (see Table H). Exchange currents were not explicitly 
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included in previous analyses of this problem. 



Table 2. The transverse electric quadrupole form factor G£;2{c|^ = 0) for the 7 + p — » A+ 
transition calculated with (i) the one-body charge density pj^^j using Siegert's theorem (first 
row), (ii) with tlie spatial current density J = J[i] + J[2] (last row) for various quark models. 
A comparison of the results in the first and last rows shows that the continuity equation is 
approximately satisfied- provided that the spatial two-body exchange currents J p] required 
by Eq.(^) are included Q. The remaining discrepancy between theory and recent experiment 
is removed by including pp] (see Eq.(|l3|)). 
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0.0203 


0.0796 
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0.0165 




0.0118 


0.0092 


0.0076 


0.0027 


0.0058 


G'£;2(J[2]) 


0.0084 


0.0114 


0.0561 


0.0127 


0.0105 


G£2(J) 


0.0202 


0.0206 


0.0637 


0.0154 


0.0163 



4 Double Spin Flip N A Quadrupole Transition 

As we have seen, a calculation using the one-body charge operator pjx] ^ind 
Siegertjla theorem yields substantially larger values for Ge2 as pointed out 
in Ref.B. We understand now why this is so. Most of the pr^] contribiition 
to Ge2 comes from the two-body currents J [2] (see Table and Ref.El for 
further explanation). This finding suggests a different interpretation of the 
deformation of the nucleon in the quark model. However, there is still a 
large discrepancy with experiment when only is taken into account. We 
have shown that a double spin flip transition based on the two-bndy charge 
operator pp] gives values for Ge2 in agreement with experiment D. This is 
explained in the following. 

The Coulomb quadrupole operator entering Eq.(|^) is: 

T^'m^^ J dl7,p(q)r„2(q), (8) 

where p = -I- p[2] is the total charge operator. The quadrupole operator 
projects onto the ^^(q), i.e., the quadrupole component in p(q). We discuss 
the quadrupole components in and pp] separately. After a multipole 
expansion of the spin-independent one-body quark charge p[i] we see that 
the Coulomb quadrupole operator is proportional to a second rank spherical 
harmonic 

T^^(p[i])(xy2(r,), (9) 
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S=l/2 S=3/2 




L=0 L=0 

Figure 3. Double spin flip due to, e.g., the two-body gluon exchange charge operator Pp^' of 
Fig.^b). The probabihty amphtude for this 5 — » S transition is proportional to oabs ~ 1, 
and thus considerably larger than the single-quark E2 transition amplitude in Eq. (U), which 
is proportional to the small D wave admixtures and b^) . A similar diagram can be drawn 
for the two-body pion exchange charge operator Ppj'' of Fig.g(c). The Coulomb quadrupole 
strength Gc2 and the related transverse electric quadrupoia strength Ge2 calculated with 
the double spin flip is in good agreement with experiment B. 



where is a single-quark position coordinate. This is the only allowed tensor 
of rank 2 that can be constructed from a one-body operator. The operator 
in Eq.(||) has nonvanishing matrix elements between the wave functions of 
Eq.(|) only for the off-diagonal 3^0,0-^3, and the diagonal D D 
transitions. 

On the other hand, the two-body gluon and pion exchange charge densities 
P[2] contain, just like the corresponding potentials, a tensor in spin space, and 
the quadrupole operator is 

T^^(p[2])«kxa]f. (10) 

Consequently, the operator in Eq. ( px| ) has a nonvanishing matrix element also 
for an 5 — > transition. We stress that unlike the single-quark operator in 



Eq.(^, the two-body quadrupole operator in Eq.(lO) does not change the 
angular momentum of the wave function. However, as a tensor of rank 2 in 
spin space it simultaneneously flips the spin of two quarks. The probability 
amplitude for this double spin flip transition is proportional to ashs ~ 1, i.e., 
two orders of magnitude larger than the iirbital angular momentum changing 
one-body transition of Eq. (^) . Morpurgolj has anticipated the important role 
of the operator in Eq.(ll0|) for the 7p — s- A+ quadrupole transition. 

The inclusion of two-body exchange currents leads to a heretofore un- 



h98: submitted to World Scientific on February 1, 2008 7 



known relation between the mean square charge radii of the N and A 

2 2 _ 2 



I: 



(11) 

Dillon and Morpurg 00 recently derived Eg. (pi]) using a rather general QCD 
parametrization method and the assumption that three-index and loop terms 
are small. They also estimate that if three-index and loop terms are included, 
deviations from Eg. (pi]) amount to 10 — 20%. 

With the help of Eg.(pl|) one obtains a parameter-independent relation 
between the neutron charge radius and the N —> A transition guadrupole 
moment 

Qp-A+ = -^{4 - ri+) = -^rl (12) 

Our parameter-independent prediction for transition guadrupole moment 
Q^p^^+ — —0.083 fm^ is in excellent agreement with the value Q^_^^+ = 
-0.086(13) fm2 extracted from the MAMI data, and Q^if^? = -0.105(16) 
fm^ extracted from the LEGS data. Furthermore, using the empirical charge 
radius of the neutron, and Siegert's theorem, we obtain using both D waves 



an 



d exchange currents Q: 



G_E2(0) = Ge2{P[1\) +Ge2{P[2\) 

= 0.017 + 0.107 = 0.124. (13) 



This has to be compared to the experimental results which lie between G^^^ 



E2 — 
in 



0.133(20) and G^'f = 0.107(17) 111. We also mention that the Z-diagrams 
Fig.^ not only explain Ge2 but also the experimental charge radius of the 
neutron u. This indicates that both the neutron charge radius and the ^ A 
guadrupole transition moment are dominated by nonvalence guark degrees of 
freedom. 



5 Double spin flip in the deuteron 

The pion exchange induced double spin-flip term also contributes to the 
guadrupole moment of the deuteron. The guadrupole moment of the deuteron 
including the pion-pair current correction is given as 

I f°° ( 1 

^^^^Tl^yo''''"^^'^ (^V0-^-2(r) 



prrY^if^r) (^2ulir) - V2uo{r)u,{r) \ul{r)^ .(14) 
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The first term in is due to the nonrelativistic one-body charge density p[i] , 
while the second term is due to the isoscalar two-body pion pair charge density 
P[2], and Yi{x) — {e^^/x){l + l/x). Evidently, the one-body contribution 
is only nonzero if there is a nonzero D-stsde wave function U2{r). On the 
other hand, even for a pure (hypothetical) 5'-wave deuteron one obtains a 
nonvanishing quadrupolc moment due to the pion-pair exchange current term. 
The Wq in the integrand is the contribution of the double spin flip term. 
Numerically, the exchange current contribution to Qd is with 0.01 fm^ rather 
small. The dominant contribution of 0.28 fm'^ comes from the one-body charge 
density involving the I?-state in the deuteron. 

Can one understand this role reversal between D states and exchange cur- 
rents when going from a single baryon to the deuteron? For a particle of mass 
m moving in a harmonic oscillator potential there is an inverse proportionality 
between the excitation energy w, the size of the system b, and the mass m, 
namely to cx l/(m6^). Because the average distance between two quarks in 
the nucleon is approximately 1.0 fm, and the average distance between the 
nucleons in the deuteron is 4 fm, one needs 2Tiuj ~ 600 MeV to lift a quark in 
an excited D state, whereas one needs only 2'hLU ~ 4 MeV to lift a nucleon into 
a D state. Given that the pion tensor potential is similar in strength in both 
systems, one can qualitatively understand why the D-state probability in the 
nucleon is so small compared to the one in the deuteron. On the other hand, 
the amount of particle-antiparticle pairs in a system is enhanced if the system 
size is decreased. Thus, the interchanged role of exchange currents and D- 
waves in the deuteron as compared to a single baryon is mainly a consequence 
of the different size of these systems. 



6 Quadrupole Moment of the A 

Tensor forces and exchange currents also lead to a nonvanishing quadrupole 
moment of the A. With configuration mixing but no exchange currents (im- 
pulse approximation) one obtains neglecting the spiall contributions and 
with typical values for the admixture coefficients a 

QT^ = b^^ (bssbDs + -^^s'.bDs^ - -O.OST&^e^. (15) 

For the quadrupole moment of the A the symmetric D state amplitude bos is 
relevant. For b = 0.61 fm one obtains Q'^^ = -0.032fm^ ga- Using V^^'^^ 
of Eq. (|l]) , Richard and Taxil, and Krivoruchenko and Giannini found = 
(2/5)r^eA, which for b — 0.71 fm coincides with Eq. (|l5|) . 
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On the other hand, with exchange currents but no configuration mixing 
we found Q 

QT = rleA. (16) 

Inserting the measured neutron charge radius into Eq.(|l6|) gives Q^'^ = 
— 0.119fm^eA- Due to the smahness of the admixture coefficients in the A 
wave function, this result remains quahtatively unchanged if the two-body 
operator is evaluated with mixed wave functions. Thus, nonvalence quark 
degrees of freedom effectively described here by quark-antiquark pair currents 
provide the dominant contribution to the quadrupole moment of the A. 



7 Deformation of the Nucleon 



This section is somewhat speculative. Because the nucleon has J — 1/2, its 
quadrupole moment in the laboratory frame is zero. This is analogous to 
a strongly deformed J = nucleus. All orientations of a deformed J = 
nucleus are equally probable. This results in a spherical charge distribution 
in the ground state and a vanishing quadrupole moment Q in the labora- 
tory. Nevertheless, one can obtain information on the intrinsic quadrupole 
moment^ Qo by measuring electromagnetic quadrupole transitions between 
the ground and excited states, or by measuring the quadrupole moment of an 
excited state with J > 1/2 ai that nucleus. 

In the collective model t3 the relation between the observable quadrupole 
moment Q and the intrinsic quadrupole moment Qg is 

3g^-J(J+l) 

(J+l)(2J + 3)^°' ^^^^ 

where J is the total spin of the nucleus, and K is the projection of J onto the 
z-axis in the body fixed frame (symmetry axis of the nucleus). The intrinsic 
quadrupole moment Qq characterizes the deformation of the charge distribu- 
tion in the ground state. A simple model for a nonspherical homogeneous 
charge distribution is that of a rotational ellipsoid with charge Z, major axis 
a along, and minor axis h perpendicular to the symmetry axis: 

07 4 
Qo = —{a^^h'')^-ZRH, (18) 

with 5^2{a- b)/{a + h) and i? = (a + 6)/2. 

"^Here, intrinsic means with respect to a body-fixed coordinate system that rotates with 
the nucleus. 
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Can one use this model to determine the sign and magnitude of the de- 
formation of the nucleon? Because the quadrupole moment of the A and the 
— + A transition quadrupole moment are dominated by gluon and pion de- 
grees of freedom, and not by single quark degrees of freedom, the collective 
model maybe valid. Assuming that the collective model is applicable, we con- 
sider the A with spin J = 3/2 as a collective rotation of the entire nucleon 
with an intrinsic angular momentum K — 1/2. One then finds from Eq.(p^) 
and Eq.([l^) Qq = 0.565 fm^ for the intrinsic quadrupole moment of the nu- 
cleon. Furthermore, from the naive nucleon model of Eq.(|l^) one obtains with 
R — yjb/irp — 1.113 fm a deformation parameter 5 ~ 0.57 and a ratio of 
major to minor semi-axes a/b « 1.79. This magnitude of deformation seems 
to be quite large. Yet, we speculate that the sign of the intrinsic quadrupole 
moment given by Eq.(|l7|) is correct. If so, the nucleon is a prolate spheroid. 
The A is clearly an oblate spheroid with an observable quadrupole moment 
given by Eq.(|6|). 

Although the nucleon wave function is not observable, one prefers to de- 
scribe the deformation of the nucleon in terms of a D state admixture in the 
quark wave function. If both S and D wave amplitudes have the same (op- 
posite) sign, one obtains a positive (negative) intrinsic quadrupole moment, 
and a prolate (oblate) deformation, respectively. For example, the sign of the 
relevant symmetric Ds wave amplitude bos in tliG ^ YftYfi function is oppo- 
site to the dominant S'-wave amplitude 65,, (see Refs. 13'BEl). This leads to a 
negative quadrupole moment of the A+ and an oblate deformation of the A, 
as can be seen from Eq. (^5|) . Because of the negative D wave amplitude in the 
nucleon wave function, one could conclude that the nucleon has an oblate de- 
formation. However, one would then miss the contribution of the nonvalence 
quark degrees of freedom. One can (via a unitary transformation) eliminate 
the two-body transition operator, as discussed by Friar for the deuteron, 
and reexpress it in terms of a modified potential, or modified wave functions. 
The matrix element is thereby left invariant. If the nucleon has an intrinsic 
quadrupole moment Qq = 0.565 fm^, the nucleon wave function should have 
a large and positive D wave amplitude, in contrast to the small and negative 
D-state admixture obtained with the original potentials. 

8 Summary 

With the help of Siegert's theorem we have shown that the largest part of 
the E2 transition strength Ge2 in — > A comes from a two-body spin 
flip due to exchange currents. The effect of one-body currents on Ge2 is 
relatively small. Exchange currents also dominate the quadrupole moment 
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of the A(1232) U,j^d are very important for the radiative E2 decays of all 
decuplet baryons llj. The quadrupole moment of the A together with the 
N —^ A transition quadruple moment can be seen as an indication for an 
intrinsic nucleon deformation. We conclude that the intrinsic deformation 
of baryons lies mainly in the nonvalence quark degrees of freedom {qq pairs, 
pions, gluons). Our prediction of the E2 transition amplitude, which is based 
on the parameter-independent relation of Eq.(p^) and Siegert's theorem, is in 
very good agreement with experiment. 

Acknowledgement: I thank H. Clement, D. Drechsel, and A. Faessler for 
stimulating discussions. 
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